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Abstract 



The effects of multiple interactions in colliding particles (e.g. in nucleus-nucleus colli- 
sions) are modeled using the light-cone dipole approach. Guided by the abelian analogue 
of multi-photon interactions in the production of a pair of charged particles, we relate the 
inclusive cross section of quark pair production with the cross sections of interaction of a 
QCD dipole with either the beam or the target. 



a,: 

1 Light-cone dipole representation 

The light cone dipole description of hadronic interactions |TJ |2| offers quite an effective phe- 
nomenology. The central quantity of this approach, the universal and flavor independent cross 
section of interaction of a colorless dipole (quark-antiquark, or glue-glue) with a target (pro- 
ton) is fitted to data, and therefore it incorporates information on all possible gluonic exchanges 
and bremsstrahlung including also nonperturbative effects. This may be considered as an al- 
ternative to the parton distribution function, with the advantage that it includes by default all 
higher order corrections and higher twist effects. This approach is especially powerful for cal- 
culating nuclear effects and diffractive processes [TJ ED EJ |3] . Since QCD dipoles are eigenstates 
of interaction, multiple interactions effects for the elastic amplitude can be included via simple 
eikonalization. 

The main difficulty of the dipole approach, unresolved so far, is modeling the distribution 
amplitude of QCD dipoles in the projectile high energy particle. This problem has only been 
solved in the lowest order of perturbative QCD for a photon projectile jHllZl, and for radiation 
of photons and gluons by a color charge jSl El UD]- However, once the multi-gluon exchange 
interactions with the target are important, one should not restrict oneself to a single parton 
density in the projectile, as is frequently done. Inclusion of higher order corrections and soft 
multiple interactions in the projectile particle remains a challenge. A model for simultaneous 
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inclusion of multiple interaction effects both in the beam and target was constructed in [11] for 
gluon radiation in nucleus-nucleus collision. Here we present another attempt to make progress 
in this direction. 

The paper is organized as follows. In Sect.|2]we study the production of a pair of charge par- 
ticles employing abelian dynamics. We start with the Born approximation (Sect. EH} and then 
include multi-photon exchanges with the target, and eventually with both colliding particles A 
and B (Sect. 12^ . 

In Sect. H3 we develop the formalism for nonabelian dynamics. The Born approximation 
is described in Sect. 13. 1| and the results are generalized to include multiple interactions in 
Sect. E21 



2 QED analogue 

2.1 Born approximation in QED 

We start with an abelian analogue for quark pair production, since this process has a simpler 
dynamics, but it contains many features of the nonabelian description. The process under 
discussion is production of a pair of particles, 1 and 2, in the collisions of two hadrons (or 
nuclei), A and B, 

A + B^A + B + 1 + 2. (1) 

For the sake of simplicity we assume that the produced particles are spinless. The cross section 
of this process reads, 

da A B = \M\ 2 dT (2) 

where M is the Lorentz-invariant amplitude, and dT is a phase space element. In the lowest 
order of perturbative expansion corresponding to the graphs shown in Fig. [IJi,b the amplitude 




Figure 1: Production of charged particles 1 and 2, in the Born approximation 
(a and b), and including multi-photon exchanges with the target or with both 
colliding particles (c and d respectively). 
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has the form, 

m = 4 a Kpa, u) 4 b) (pb, to ) v(^y, P2 ) ( (3) 



where 

qA(B) = Pa(b) - p'a(b) ! 

Qa(B) = -QA(B) ~ 9A(B) ± + lA(B) min 5 

iA(B) mm = -j^- ; (4) 

7a(b) 

Pa(b), Pa(b)i Pi an d £>2 are the 4-momenta of particles in the initial and final states, and of 
the produced particles respectively; ja(b) are the Lorentz factors of the colliding particles; qA(B) 
are the 4-momenta of the virtual photons emitted by the particles A(B). The electromagnetic 
currents in Eq. (JHJ) have the form, 

4 A) (PA, q A ) = V^^Z A F^\Q 2 A )(2p A - q A )^ ; (5) 
4 A) (PB, q B ) = V^ZZ B F^ B \Q 2 B )(2p B - q B ) v , (6) 

where Z A (b) are the charges of A and B, and F A<yB \Q\^) are their form factors, respectively 
The tensor in Q reads, 



D 1 1 (2pi - q A )^{2p 2 ~ q B )v + D 2 1 ( 2 P2 - <?a) m (2pi - q B ) v - ; (7) 

D 1 = m 2 - ( Pl - q A f = m 2 - (p 2 - q B f ; (8) 

D 2 = m 2 - (p2 - g A ) 2 = m 2 - (p x - g B ) 2 ; (9) 

Notice that, 

(?A)/i*i» = V (to)* = • ( 10 ) 
The phase space factor in Eq. (J2J) is 



d 3 p' A d 3 p' B d 3 p 1 d 3 p 2 (4) , , 

^ = 64 7 (2VT)* £^ ei 6, 5 iPA ~ ~ ~ P 1 ~ P») 

64J(2vr) 8 eie 2 



64/(2 7 r) 8 J B^e 1 e 2 



where 



( P aPb) 2 - M 2 M 2 . (12) 



Here Ma(b) and £?a(b) are the colliding hadrons or nuclei masses and energies; m and e^) are 
the masses and energies of the produced particles. 

To simulate QCD, we introduce a charge screening effect in what follows, i.e. consider the 
colliding particles as neutral dipoles with the screening potential, 

K(r) = ^e" Ar . (13) 
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This can be also simulated by an effective photon mass A, replacing in (JHJ), Qa(b) ^ Qa(b) + ^ 2 - 
Now, let us consider the collision of ultra-relativistic particles ('Ja(b) ^ 1) i n the cm. frame. 
Neglecting corrections of order 0(7^^) we have from © - (H2J), 



a 



dq\dq A dq\dq B 
+ 



d 3 pid 3 p2 A/ \ B 1 \\TTI \|2 

{QA)cr {q B ) \U(qA,qB,Pi,P2)\ 



£1 e 2 



x S(Qa) 6 (1b) <K<?a + Qb ~Pi 



Pt) S(q~, 



q B ~Pi ~ P2 



(14) 



where pf^ = (pi(2))o±(pi(2))« are the light-cone momenta of particles 1(2) (the z-axis is chosen 
along the momenta of A, B); a^q^.) and cr®(q B ) are the differential cross sections of elastic 
scattering of particles 1, 2 on hadrons (nuclei) A or B respectively, 



<4(u) = 

Then, we can represent 



d 2 a[l{2) + A ^ 1(2) + A] 
d 2 q A± 

d 2 a[l(2) + B -> 1(2) + 5] 



d 2 q 



d 3 pi 



(2) 



dpi 



(2) 



d 2 p m± 

£ 1(2) Pl(2) 

+ - _ 2 I —*2 
Pl(2)Pl(2) - m +Pl(2) ± 



d 2 pi(2)„ 



2Z A a em F A (Q 2 A y'' 

+ A2 

~2Z B a em F B (Q 2 B 

+ A2 

Pl(2) 



(15) 
(16) 



(17) 
(18) 



It is convenient to introduce the fractions of light-cone momenta of the colliding virtual 
photons carried by the produced particle 1, 



Pi 

UA = — r 



a B = — 



Pi_ 



which are connected by the relation, 

<*A{B) 



(1 - atB(A))(mi)i + a B(A ))(m 2 ) 
The amplitude U(q A ,q B ,pi,p 2 ) in (H2)? 



(19) 



(20) 



U(q A ,qB,Pi,P2) 



Pi P2 , P2 Pi 
Di 



Dn 



(21) 



is function of three 4-momenta (since q A + Qb — Pi + ^2)- Therefore, we can choose as inde- 
pendent variables the three transverse momenta and one of the light-cone fractions a A (or a B ). 
Selecting q A± , q B± , p 1± and a^, we can represent U as, 



U(qA ± ,qB ± ,Pi ± ,at A ) = ${Pi x ~ a A q A± ; q A± , a A ) - $(pi ± - a j4 ^ i - qB ± ] U ± , a A ) 



(22) 
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where 



®(p±;q±,a) = ^(p±;q±,a) + ^ L (p ± ;q ± ,a) ; (23) 
P± + e 2 {q ± ,a) 

,^ ^ a(l — a)(l — 2a) g? /ri _ N 

p± + e 2 {q±,a) 

e 2 (q±,a) = m 2 + a(l — a)q]_ . (26) 
Thus, the cross section of pair production can be represented as, 

a(A + B ^ A + B + 1 + 2) = J dq\ d 2 q A± n A (q A ) a( 7 * A +B -> 1 + 2 + B) , (27) 

where 

a. 



n ^) = FoZvii I „+ = I ^2 I ^ ) — > (2* 



(27r)2a em gi WI + A V OX 

is the density of equivalent photons |T21 H3| in the projectile A. 

The virtual photoproduction cross section a{^* A + B — > 1 + 2 + 5) in ()27j) can be expressed 
in terms of the dipole formalism as, 

a( lA + B^ 1 + 2 + B) = J d 2 rda A [\^ T A (f,a A )\ 2 + \^(r,a A )\ 2 ] a B (r) . (29) 

Here \I/ T ' L are the light-cone wave functions of transversely or longitudinally polarized photons 
with 4-momentum q A , 



a A ) = J " em . F ,. o / d 2 p T ^(pl; gk x , a A ) e*^ . (30) 



a A (l - a A ) \q Al \(2ir) 

The cross section of interaction of the dipole of particles 1 — 2 with B has the standard form 

in, 

a B {f) = 2 J d 2 q B± a*(q B ) (l - e^^) . (31) 

The resulting representation, Eq. (J29)) . which treats the production of particles 1 — 2 as 
photoproduction by a virtual photon representing the electromagnetic field of the projectile A 
interacting with the target B, looks asymmetric relative to the replacement A ^ B. This is, 
however, an artifact of our choice of a A as a variable. If our choice were a B , the same cross 
section would look as a result of interaction of a photon 7^ with the target A. Thus, the choice 
of an independent variable, a A or a B , leads to a breaking of the symmetry, A ^ B. 

2.2 Multi-photon exchanges 

So far our considerations were restricted to the Born (one photon) approximation. It turns out, 
however, that the relation Eq. (J29|) is also correct if the dipole 1 — 2 interacts with the target 
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B via multiple-photon exchanges as is illustrated in Fig. ^p. This is particularly important if 
B is a nucleus. Indeed, in this case the dipole cross section reads, 



a B {r) = 2 J d 2 b{l -exp[zAx B (&,r)]} , (32) 



where the phase shift 



Axb(M = XB (b + ) - XB {b-) ; (33) 

6+ = b + (l-a A )r, 
= b — a a? . 

At first glance the cross section Eq. ()32j) depends also on a^. However, the change of integration 
variable, b =>- b + (a a — l/2)r eliminates the a a dependence. 

Another possible representation for the dipole cross section has the form, 

a B (r) = 2 J d 2 q ± ag(q ± ) (l - e^) , (35) 

where cr^(gl) is the differential cross section of elastic scattering of one of the particles, 1 or 
2, with the target B, calculated in the eikonal (Glauber) approximation, 



°gi(Q±) = ( 36 ) 



2 



/±(£l) = ± I d 2 be^ b l- e ±^B(6) . (37) 
2n J L - 1 

The signs ± correspond to the opposite charges of particles 1, 2. 

Thus, we conclude that the effect of all multi-photon exchanges with the target B (restricted 
to only one photon exchange with the projectile A) is equivalent to the replacement of the Born 
cross section cr B (q B ) by the Glauber one, cr^(g B ). 

The same cross section of pair production calculated in a single photon approximation for 
A, but multi-photon with B, can be represented differently if one chooses a B as a variable, 

a( 7 *a +B -> 1 + 2 + B) = J dq B _ d 2 q B± n B {q B ) o( 1b + A^l + 2 + A). (38) 

Here the process j B + A — > 1 + 2 + A is calculated in one-photon approximation, but the density 
function of "equivalent photons" is different from 



B 



{qb) Qb i 



n B (q B ) = t^ttP , (39) 
(27r) 2 a em q B 

with the replacement of single- to multi-photon exchange, cr B (q B ) =>- ctqAQb)- 

It is natural to assume that the inclusion of multi-photon exchanges between the produced 
quark pair and both colliding nuclei (Fig. [TJI) can be done by replacing the " single-photon" 
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quantity a{^* B + Z A — > 1 + 2 + Z^) in (J3S|) by a multi-photon one. To do that we should replace 
the Born approximation for the cross section a A (r) in (}2"9"|) (with the interchange A <^ B), 



a A (r) = 2 J d 2 q A a£(q A ) (l - e^) 



by the following Glauber form, 



a A (r) = 2 d 2 b 



(40) 



(41) 



The phase shifts Ax A (b, r) are defined on analogy to Axb(6, f) in (J3HJ)-(J3IJ)- 



3 Quark production 

3.1 Born approximation in QCD 

The Feynman graphs corresponding to qq pair production in a collisions of two hadrons A 
(beam) and B (target) in lowest order in a s (double one-gluon approximation) are shown in 
Fig. El We assume that only one quark (antiquark) is detected, with transverse momentum 



A 




a b c 

Figure 2: One gluon approximation to the central production of a qq pair. 



and rapidity y, while the accompanying antiquark (quark) is not observed, i.e. its momentum 
is integrated out. Then the cross section corresponding to the graphs in Fig. El has the form, 



da (A B — >• qX) Akol s f dx\ 



d 2 pr dy 



Id q 1 d q 2 



q 2 J r A (x 1 ,q 1 ) 

mini 



— (7 $1 

16 1 1 



91$; 



^B(x 2 ,q 2 ) 

2 \2 
mini 



(?2 + ?2 



(42) 

Here pr and y are the transverse momentum and rapidity of the produced quark (or antiquark); 
x q is the fraction of the plus component of the momentum of the hadron h\ taken by the quark, 
which is related to the rapidity interval Ay = ln(l/a; g ) between the hadron and the quark; 
and qip are the transverse momenta of the gluons radiated by the hadrons h\ )2 with light-cone 
fractional momenta X\ t2 respectively. While we integrate over x%, the value of x 2 is defined by 
the kinematics, 



x 2 



1 

XiS 



a 



+ 



m l + (pt -qi- &) 
1 — a 



(43) 
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where a = x q jx\. The function <3>! in (J42)) corresponds to the sum of diagrams Fig. |2t,b, while 
$2 corresponds to the difference of the amplitudes Fig. |2^,b plus the graph in Fig. |2b. These 
functions are expressed in terms of the usual LC wave functions (fcy, a, q 2 ) of the qq Fock 

state in a gluon, where kx is the relative transverse momentum of the qq, 
where 

e 2 = a(l - ct)Q 2 + m 2 ; 

Ql = Ql+ Qlmin (Qmin ~ ^-QCd) ] 

O = ma ■ e + (1 — 2a) (a ■ n)(p ■ e) + i(p x n) ■ e ; (45) 

and the gluon polarization vector is related to its transverse momentum, e = q\jq\- Then, we 
have 

$i = ^qqipr ~ aqi, a, q\) - ^%{p T - aqi - ql, a, q\) ] 

*2 = ^%(PT ~ aqi, a, q\) + (p T - aq\ - <f 2 , a, g 2 ) - 2 (p T - aq\ - aq 2 , a, q\). (46) 



3.2 Multiple interactions 



In the light-cone approach employed in the rest frame of the target (the bottom hadron in Fig. [2J) 
the process depicted in Fig.|2]looks like the interaction of a qq fluctuation of the projectile gluon 
with the target. Although the interaction is mediated by one gluon exchange, one can make 
it more realistic using the phenomeno logical dipole cross section a^ q (r,x) of interaction of a qq 
dipole of transverse separation r with a hadron h at energy s ~ (xr 2 )~ l . This cross section 
fitted to data incorporates the unknown dynamics of soft multi-gluon exchanges and radiation. 
Applying to ()42|) a Fourier transformation we get, 



da(AB -> qX) 
d 2 px dy 



1 r dxA \d\ x d\^ TA{x ^ 



57T 



X\ 



x 



d 2 r x d 2 r 2 exp i(p T - aq 2 )(n - f 2 ) *f (r 2 , a, q\) *£(fi, a, g 2 ) £ B (r*i, f 2 , x 2 , a), (47) 



where 



E (n,r 2 ,x 2 , a) 



16 



af 2 ,x 2 ) + <T^(f 2 



ar\,x 2 ) + crf q (afi,x 2 ) + crg a [ar 2 . r 2 ] 



vf q {ri,x 2 ) 



Vq q {r 2 ,x 2j 



r 2 ,x 2 



■qq\ 



ar 2 ,x 2j 



}. (4* 



At r*i = r 2 this becomes the familiar combination cr 3 (r, a, x) = |{((T^(ar, x) + a q 1 q [(l — a)f, x}} — 
|cr^[(2a — l)r, x] which is the dipole cross section for a three-body system qqG) interacting 
with a hadron h. In particular, it enters the total cross section of qq pair production by a gluon 

us Gang. 



8 



The dipole cross section vanishes at small qq separations, cr£ (r, x) rT ^ = r 2 G h (x, 1/r 2 ) Tr 2 a s /3, 
where G h (x, 1/r 2 ) = xg h (x, 1/r 2 ) is the gluon distribution function in hadron h. In this limit the 
dipole cross section corresponds to one gluon exchange (in the inelastic amplitude) and Eq. (}4*2*|) 
is recovered. At the same time, it is usually assumed that at large separations a^ q (rT,x) sat- 
urates at some constant value ctq. This may be motivated by either saturation of the gluon 
density, or shortness of the gluon interaction radius. 

Let us introduce a function 



uJ h {r,x) -<4(r» , 

which has the properties u h (f, x) r ^ — ► <7q(x) and u(f, x) r „ 
transform, 

u h (q,x) 



(2tt)' 



d z ruj h {r,x) e iq - r 



(49) 

0. Therefore its Fourier 
(50) 



is defined for any q. Then, taking into account that / d 2 qu} h (q,x) e %q ' T = <Jq(x), we can 
represent the function £ B (ri, f 2 , a), Eq. (jlHj) . in the form, 



S (f 1 ,f 2 ,x 2 a) 



+ 



9 
16 



1 + e tq - ri - 2 e 



iaq-ri 



1 + e iq ' r2 - 2 e 



iaq-T2 



(51) 



Using this expression we can rewrite Eq. (|37jl as 



da(AB^qX) 1 



d 2 pT dy 



a 



- / dx\ — Id qi 



Xi 



2 q 2 J z A {x 1 ,q l ) B 
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.x 2 ) 



7 ,2 9 ._ |2 
— $1 H $2 

16 11 16 1 1 



(52) 



Finally, using the relation ui(q,x) = ao(x)d(q) — (jq q (q,x) we arrive at, 



da(AB^qX) 1 
d 2 px dy 2 



dx\ 



a 



d 2 qi 



glF A {xi,qi) / 7_ 

(Ql+Qlmin) 2 V 16 



i$ir + 



9 

19 



I $9 



(53) 



Comparing this expression with Eq. f!42D we conclude that the transformations done above are 
equivalent to the replacement 



(&Fb{x<iAi) 3 B 



(54) 



in Eq. (J42)) . This observation leads to the natural assumption that the same procedure should 
be performed with the contribution to Eq. (jUJ) of the upper part of the graphs in Fig. El namely, 



q 2 J r A (x 1 ,q 1 ) 

mini 

Switching back to coordinate representation, 



47ra;. 



(27T) 



(55) 



(56) 



we eventually arrive at the cross section in a form which includes multi-gluon exchange, 



da(AB -> qX) 
d 2 p T dy 
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(47r) 3 a s 



dxi 



a 



d 2 q / d 2 pd 2 r 1 d 2 r 2 exp i(pr — aq) ■ (fi — r 2 ) — iq - p 



(57) 



x Vjc4(p, xi) tf^fa, «, q 2 ) *g(fi, a, g 2 ) S B (f 1; f 2 , x 2 , a) . 

This is the central result of this paper. It looks asymmetric, while the graphs in Fig. El are 
symmetric relative to beam-target interchange. This expression, however, has been derived in 
the rest frame of the target. In the beam rest frame one should just interchange A <^> B and 

X\ <^> X 2 . 

The total yield of quarks integrated over transverse momentum has the simple form, 



J dxx^ J d 2 q J d 2 pd 2 re~ i '* PV 2 p (Tf q (p,x 1 ) ^?(f,a,g 2 



.B, 



da(AB -> qX) _ 6 

dy (47r) 3 a s 

(58) 

Thus the production cross section is expressed in terms of the dipole cross section either on 
one (A), or another (B) colliding particles or nuclei. Unfortunately, it has an asymmetric form 
which is related to our choice of a^. Switching to «g we will get an equivalent cross section, 
but having a different form. 



4 Summary 

Guided by the abelian analogue of particle production in QED, we suggested an approach incor- 
porating multiple interactions both in the beam and target, within the dipole formalism. The 
proposed procedure of replacing the unintegrated gluon density in the target by a combination 
of cross sections of dipole-target interaction is assumed to be valid also for the beam. The main 
result, Eq. (JH7|) . still needs needs to be tested through numerical calculations and comparison 
with data. This expression looks asymmetric relative to the dipole cross sections of interacting 
with the colliding particles or nuclei, A and B. However, effectively it is symmetric provided 
that a proper replacement of variables is done. It is possible to rewrite it in an explicitly 
symmetric form, however, in this case the light-cone distribution function ^>q q in (|H7j) should 
be replaced by a more complicated function. Further development of this formalism will be 
published elsewhere. 
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